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Multiple parton scattering and induced parton energy loss in deeply inelastic scat- 
tering (DIS) off heavy nuclei is studied within the framework of generalized factoriza- 
■ tion in perturbative QCD with a complete calculation beyond the helicity amplitude 

(or soft bremsstrahlung) approximation. Such a calculation gives rise to new cor- 
f"^ ' rections to the modified quark fragmentation functions. The effective parton energy 

CNj , loss is found to be reduced by a factor of 5/6 from the result of helicity amplitude 

(-H ■ approximation. 
C3 
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7-H ' I. INTRODUCTION 

> : 
^ ■ 

^\ . Suppression of jet production or jet quenching in high-energy nuclear collisions has been proposed 

as a good probe of the hot and dense medium [1,2] that is produced during the violent collisions. 
The quenching of an energetic parton is caused by multiple scattering and induced parton energy 
loss during its propagation through the hot QCD medium. It suppresses the final leading hadron 
distribution giving rise to modified fragmentation functions and the final hadron spectra [3,4] . Recent 
theoretical estimates [5-9] all show that the effective parton energy loss is proportional to the 
gluon density of the medium. Therefore measurements of the parton energy loss will enable one to 
extract the initial gluon density of the produced hot medium. Strong suppression of high transverse 
. momentum hadron spectra is indeed observed by experiments [10,11] at the Rclativistic Heavy-Ion 

•rH , 

Collider (RHIC) at the Brookhaven National Laboratory (BNL), indicating large parton energy 
loss in a medium with large initial gluon density. However, one cannot unambiguously extract the 
initial gluon density from the experiments of heavy-ion collisions alone because of the theoretical 
uncertainty in relating the parton energy loss to the initial gluon density. For this purpose, one has 
to rely on other complimentary experimental measurements such as parton energy loss in deeply 
inelastic scattering (DIS) of nuclear targets. One can then at least extract the initial gluon density 
in heavy-ion collisions relative to that in a cold nucleus [12]. 

Modified quark fragmentation function inside a nucleus in DIS and the effective parton energy 
loss has been derived recently by Guo and Wang [13]. Generalized factorization of twist-four pro- 
cesses [14] was applied to the inclusive process of jet fragmentation in DIS in order to derive the 
modified fragmentation functions. Taking into account of gluon bremsstrahlung induced by multiple 
parton scattering and the Landau-Pomeranchuck-Migdal (LPM) interference effect, one finds that 
the leading twist-four contributions to the modified fragmentation function and the effective parton 
energy loss depend quadratically on the nuclear size Ra- They also depend linearly on the effective 
gluon distribution in nuclei. One can also extend the study to parton propagation inside a hot QCD 
medium reproducing earlier results [12]. This allows one to relate parton energy loss in both hot 
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and cold nuclear medium. 

There are all together 23 cut-diagrams that contribute to the leading twist-four corrections to the 
quark fragmentation function in eA DIS. For simplification of the calculation, the helicity amplitude 
approximation was used in Ref. [13] in the limit of soft gluon radiation z g = 1 — z — > where z g 
is the momentum fraction carried by the radiated gluon and z the fraction carried by the leading 
quark. Such an approximation enables one to simplify the calculation of the radiation amplitudes. 
The final results are obtained by squaring the sum of all possible amplitudes, giving rise not only 
to the contributions of double scattering but also various interferences. In this approximation, the 
amplitudes of initial and final state radiation are the same except the opposite signs and different 
color matrices. Because of the different color matrices in the initial and final state radiation, there is 
no complete cancellation of the radiation amplitudes. In addition, there is also gluon radiation from 
the exchanged gluon via triplc-gluon coupling. These non-Abelian features of QCD radiation lead 
to a finite gluon spectra even in the helicity amplitude approximation. However, under the same 
approximation, the photon spectra from QED bremsstrahlung would be zero because of almost 
complete cancellation between initial and final state radiation. One therefore has to go beyond 
the helicity approximation. In this paper, we will study the correction to the gluon radiation 
spectra when we go beyond the helicity amplitude approximation and its effect in the modified 
quark fragmentation function. We will also compute the effective quark energy loss and compare to 
the result in the helicity amplitude approximation. 



II. GENERALIZED FACTORIZATION 



In order to study the quark fragmentation in eA DIS, we consider the following semi- inclusive 
processes, e{L\) + A(p) — ► e(L 2 ) + h(£h) + X, where L\ and L 2 are the four momenta of the 
incoming and the outgoing leptons, and lh is the observed hadron momentum. The differential cross 
section for the semi-inclusive process can be expressed as 

wherep= [p+,0,Oj_] is the momentum per nucleon in the nucleus, q = L 2 — L\ = [— Q 2 /2q~, q~ , 0j_] 
the momentum transfer, s = (p + Li) 2 and o;em is the electromagnetic (EM) coupling constant. The 
leptonic tensor is given by L^ v = l/2Tr(7 • £17^7 • L2lv) while the semi-inclusive hadronic tensor 
is defined as, 



dH h 2 

x2Tr5 4 (q+p- Px -£ h ) (2) 



x 
nit 



where ^2 X runs over all possible final states and J M = J2 q Zq^qlix^q is the hadronic EM current. 

In the parton model with collinear factorization approximation, the leading-twist contribution 
to the semi-inclusive cross section can be factorized into a product of parton distributions, parton 
fragmentation functions and the partonic cross section. Including all leading log radiative corrections, 
the lowest order contribution (0(a°)) from a single hard 7* + q scattering can be written as 
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dW s f 

= E e 'y dxf^x, f , 2 I )HjZ ) (x 1 p,q)D q ^ h (z h ^ 2 ) ; (3) 

1 27T 

H$(x,p, q) = - Tr( 7 • ■ (q + xp) lv ) ^—8(x ~ x B ) , (4) 

where the momentum fraction carried by the hadron is defined as zj- L = 17 /q~ and xb — Q 2 /2p + q~ 
is the Bjorken variable, /if and p? are the factorization scales for the initial quark distributions 
f q (x, p, 2 ) in a nucleus and the fragmentation functions D q ^i l (zf l , /i 2 ), respectively. The renormalized 
quark fragmentation function D q ^h{zh, p 2 ) satisfies the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi 
(DGLAP) QCD evolution equations [16]. 

In a nuclear medium, the propagating quark in DIS will experience additional scatterings with 
other partons from the nucleus. The rescatterings may induce additional gluon radiation and cause 
the leading quark to lose energy. Such induced gluon radiations will effectively give rise to additional 
terms in the evolution equation leading to the modification of the fragmentation functions in a 
medium. These are so-called higher-twist corrections since they involve higher-twist parton matrix 
elements and are power-suppressed. We will consider those contributions that involve two-parton 
correlations from two different nucleons inside the nucleus. They are proportional to the size of the 
nucleus [17] and thus are enhanced by a nuclear factor A 1 / 3 as compared to two-parton correlations 
in a nucleon. Like in previous studies [13], we will neglect those contributions that are not enhanced 
by the nuclear medium. 




L C R 

FIG. 1. A typical diagram for quark-gluon re-scattering processes with three possible cuts, central(C), 
left(L) and right(R). 



We will employ the generalized factorization of multiple scattering processes [14]. In this ap- 
proximation, the double scattering contribution to radiative correction from processes like the one 
illustrated in Fig. 1 can be written in the following form, 



dz h 



a0 



d 2 — D 



1 



2 dk%dl4 " " 



(5) 



J k T =0 



after collinear expansion of the hard partonic cross section with respect to the transverse momentum 
of the initial partons, where H „„ (y~ ,kx,p, q, z) is the Fourier transform of the partonic hard 

part H^ v (x, xi, X2, kr,p, q, z) in momentum space, 
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H^(y-,y^,y^,k T ,p,q,z) = J dx ^ ^ e ix lP +y-+ix 2P +yr+i(x- Xl - X2 ) P +y; 

xH® v (x,xi,x 2 ,k T ,p,q,z) , (6) 

and kx is the relative transverse momentum carried by the second parton in the double scattering. 
This is the leading term in the collinear expansion that contributes to the double scattering process. 
The first term in the collinear expansion gives the eikonal contribution to the leading- twist results, 
making the matrix element in the single scattering process gauge invariant, while the second (or 
linear) term vanishes for unpolarized initial and final states after integration over fcy- 

The hard part of the partonic scattering for each diagram, H^ u (x, xi,X2, kr,p, q, z), always con- 
tains two (5-functions from the on-shell conditions of the two cut-propagators. These <5-functions, 
together with the contour integrations which contain different sets of poles in the un-cut propaga- 
tors, will determine the values of the momentum fractions x,x\, and X2 [13]. The phase factors 
in H^ l/ (y^ , 1 y^ , kx,p, q, z) [Eq. (6)] can then be factored out, which will be combined with the 
partonic fields in Eq. (5) to form twist-four partonic matrix elements or two-parton correlations. The 
double scattering corrections in Eq. (5) can then be factorized into the product of fragmentation 
functions, twist-four partonic matrix elements and the partonic hard scattering cross section. 



III. BEYOND HELICITY AMPLITUDE APPROXIMATION 



To simplify the calculation of various cut-diagrams of double scattering and illustrate the underly- 
ing physical processes, helicity amplitude approximation was used in Ref. [13]. In this approximation, 
one neglects the transverse recoil induced by the scattering and consider only the part of the ampli- 
tudes in which quarks' helicity is unchanged in the scattering. The final results will agree with the 
complete calculation in the limit of soft radiation. 

Take photon bremsstrahlung for example. A complete calculation of photon radiation induced by 
a single scattering with transverse momentum transfer kr gives a spectra 

dN a 
dl\dz ~ 2^ 

Here we denote the momentum of the photon (or gluon in QCD) to be I which carries 1 — z momen- 
tum fraction of the struck quark. Under helicity amplitude approximation, the splitting function 
will become 2/(1 — z) and furthermore the term (1 — z)kx in the final state radiation amplitude 
will be neglected. The interference between initial and final state radiation will effectively reduce 
the photon radiation spectrum to zero, apparently not a precise approximation. In QCD, the corre- 
sponding gluon radiation has similar amplitudes, except additional color matrices. Because the gluon 
exchange in the scattering transfers color, the emitted gluon in the initial and final state radiation 
can carry different colors. In this case, there is no complete destructive interference between initial 
and final state radiation as in QED. The helicity amplitude approximation is, therefore, a better 
approximation in QCD than in QED. However, it still neglects the corrections which contribute the 
most in photon radiation in QED. This finite correction is what we will study in this paper. 



£ T e T + (i - z )k T 

4 ~ (It + (1 - z)k T ) 2 



1 + z 2 



(7) 
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We first consider the contribution from Fig. 1 in detail and will list the results of other diagrams 
afterwards. Using the conventional Feynman rule, one can write down the hard partonic part of the 
central cut-diagram of Fig. 1 [13], 



x ^Tr 
2 



) 4 

/--v, /.V7/,,r ,. (I 2 ) 6(1 - z - l —) . (8) 



~ _ Cf_ 4 7 ■ (q + xip) 7 ■ (g + xip - t) 

ap ~ 2N C 9 (q + x lP ) 2 - le la (q + x lP - If - ie lal ' lq lp 
v e aPft) l-(q + xp-£) ~f-(q + xp) 2 



(9) 



where e a/3 (£) is the polarization tensor of a gluon propagator in an axial gauge, n ■ A = with 
n = [1, 0~, 0j_], and £, £ q = q + (xi + x 2 )p + k T — £ are the 4-momenta carried by the gluon and 
the final quark, respectively, z = £~ /q~ is the fraction of longitudinal momentum (the large minus 
component) carried by the final quark. 

To simplify the calculation, we also apply the collinear approximation to complete the trace of 
the product of 7-matrices, 

^v»7-^ii[7yv] • ( 10 ) 

After carrying out momentum integration in x, xi, X2 and ^ ± with the help of contour integration 
and (^-functions, the partonic hard part can be factorized into the production of 7-quark scattering 
matrix H$ (x,p,q) [Eq. (4)] and the quark-gluon rescattering part H° , 

H^uiv ~,Vi,V2 , k T,P,q,z) = J dxH$(x,p,q) H° \y~ , ,y^ ,k T ,x,p,q,z) . (11) 

Contributions from all the diagrams have this factorized from. Therefore, we will only list the 
rescattering part H° for different diagrams in the following. For the central-cut diagram in Fig. 1 
it reads [13], 

— d _ _ _ f d£\ a s 1 + z 2 

H C (F ig .i)(y ,Vi ,Vi ,hr,x,p,q,z) = J ^Cf— 

x ~^- 1 Ic(Fig.i)(y~ ,Vi ,V2 JT,k T ,x,p,q, z) , (12) 

Ic(F ig .i)(y-, Vi, V2 At, kr, x,p, q, z) = e^ x+x ^ + y~ +ix - p+ ^ -^9(-y^)e(y- - y^) 

x (i- e - iXLP+y ')(l-e- iXLP+( ' v ~- y * ) ) . (13) 

Here, the fractional momentum is defined as 

£ 2 T _k%- 2k T ■ £ T 

XL ~ 2p+q-z(l - Z) ' XD - 2p+q-z ' [ > 

and x = xb = Q 2 /2p + q~ is the Bjorken variable. 
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The above contribution resembles the cross section of dipole scattering and contains essentially 
four terms. The first diagonal term corresponds to the so-called hard-soft process where the gluon 
radiation is induced by the hard scattering between the virtual photon and an initial quark with 
momentum fraction x. The quark is knocked off-shell by the virtual photon and becomes on-shell 
again after radiating a gluon. Afterwards the on-shell quark (or the radiated gluon) will have a 
secondary scattering with another soft gluon from the nucleus. The second diagonal term is due to 
the so-called double hard process where the quark is on-shell after the first hard scattering with the 
virtual photon. The gluon radiation is then induced by the scattering of the quark with another gluon 
that carries finite momentum fraction xl + Xd- The other two off-diagonal terms are interferences 
between hard-soft and double hard processes. In the limit of collinear radiation (xl — * 0) or when 
the formation time of the gluon radiation, t/ = 1/xlp + , is much larger than the nuclear size, the 
two processes have destructive interference, leading to the LPM interference effect. 

One can similarly obtain the rescattering part H° of other central-cut diagrams (a-d) in Fig. 2: 

di\ a. 1 + z 2 

— — (-j A 

i T _ k T ) 2 2lT 1 - z 

2na s - . 
X — Jc(a)(y ,2/l ,2/2 ,ZT,k T ,X,p,q,Z) , 

7c(a)(y-, 2/i~, y^Jr, hr, x,p, q, z) = e^ +x ^ + y~ +ix ^ + ^-y^e(-y^e(y--y^) 

x J e ^r>P + y^/(l-2) _ e -ixLP + y? I 

x [ e ix DP + (y~ -Vi)/(i-z) _ e -ix L p + (v~-yT)^ ^ 
D _ _ f di\ a s 1 + z 2 

H C ( b )(y ,y 1 ,y 2 ,kT,x, P ,q,z)= / — ; — -__c F - — - 

J (£ T - (1 - z)k T ) 2 27T 1 - z 

x ~Jf ± Ic(b){y~ ,2/i~ ,2/2 ~jT,k T ,x,p,q,z) , 
I C{b )(y-, Vi, 2/ 2 ~ ,*t, hr, x,p, q, z) = e ^ +x ^ + y~ +ix ^ + ^-y^e{-y^)6{y- - y^) 

x e -ixLP + (y~-Vi ) e -ixLP + y2 ) Mgj 

■ttD . _ _ _ , . f 2 i^r k T ) ■ [f T - (1 - z)k T ) a s C A l + z 2 



H c(a)(y ,2/i ,2/2 ,kr,x,p,q,z) = J 



H c(c)(y ,2/i ,V2 ,kr,x,p,q, z) = J dl 2 T 



(fT-k T ) 2 {f T -(l-z)k T ) 2 2^ 2 1-z 



x -jj—Ic(c)(y ,2/1 ,2/2 ,ZT,kT,x,p,q,z) , 
Ic(c) (y~ , Vi , V2 , *t, hr,x, p, q,z) = e *(*+*i> + ir-H* DP + for - Ba " ) 6{ _ y - )6{y - _ y -) 

x & -iXLP + V2 

x j e ^DP + (y _ -ar)/(i-z) _ e -ix L p + (y~ -y±)^ ; 
^ .._ .._ t , _ /■ ^ _(£ - fcr) • (£ - (1 - z)k T ) a s C A l + z 2 



H c(d)(y ,2/i ,2/ 2 ,kr,x,p,q,z) = J < 



(4r - &t) 2 (^t - (1 - z)fc T ) 2 27r 2 1 - z 



x -jj-Ic(d)(y ,2/i ,2/ 2 jT,kr,x,p,q,z) , 
lew (»", Wi" , 2/2", 4r, fer, g, z) - e^ x+x ^ + y~ +ix ^ + ^-^ 6{-y^)6{y- - y^) 
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To complete the calculation we also have to consider the asymmetrical-cut diagrams (left cut and 
right cut) that represent interferences between single and triple scatterings. They can be obtained 
with similar procedures. We list the rescattering part H° of all those asymmetrical-cut diagrams 
in the Appendix. 

To obtain the double scattering contribution to the semi-inclusive processes of hadron production 
in Eq. (5), one will then have to calculate the second derivatives of the rescattering part H° . 

After a closer examination of these rescattering parts, one can find that all contributions from the 
asymmetrical-cut diagrams have the form as 

It ■ (It ~ f(z)h 



(19) 



i\(i T ~ f(z)k T y 

where f(z) = 0, 1, 1 — z, z is only a function of z, X is the longitudinal momentum fraction and Y~ 
the spatial coordinates. One can prove that the second derivative of the above expression vanishes 
at Ut = 0, 

f T ■ (f T - f(z)k T ) 



e T (i T - f{ z )k T f 



0. 



(20) 



Therefore, all contributions from the asymmetrical-cut (right-cut and left-cut) diagrams will vanish 
after we make the second partial derivative with respect to kx when we keep only the leading terms 
up to 0{x B /Q 2 £ T ), 

tD 



V 2 H 

v Ut asym 



\k T =o = + 0{x B /Q 2 i T ) 



(21) 



In the same way we find that some of the central-cut diagrams will not contribute to the final results, 
cither. In fact, after making the second partial derivative with respect to k? only four central-cut 
diagrams shown in Fig. 2 will contribute to the final result. 
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(c) (d) 
FIG. 2. Four central-cut Diagrams that contribute to the final results. 



Including only those contributions that does not vanish after second derivative with respect to 
hr, we have 



-(1 — e - ix LP + V2 )(l — e - ix LP + (y~ -J/^)) 



Vi) 



+ 



4C F (l - zf j XL p+(y--y 7 ) 



e -ix L p^y 2 



_l_ 2Ca(1__z) e ~ix L p + y 2 ^ _ e -ixLP + (y -y 1 )-) 
2C A (l-z) . 

+ ^4 6 



■ix L p + (y -y 1 ) 



(1 



") + o(x B /<3 2 4) 



(22) 



The first term at the right-hand side in Eq. (22) comes from the contribution of #c(a) which is the 
main contribution in the previous calculation [13] with helicity amplitude approximation. It contains 
hard-soft, double hard processes and their interferences. The other three terms come from diagram 
(b),(c),(d) of Fig. 2 respectively. They constitute corrections to the first term in powers of 1 — z. 
The second term that are proportional to (1 — z) 2 is from the final state radiation from the quark 
in the double hard process in Fig. 2-(b). This term is the only contribution to the finite photon 
spectra in the corresponding QED bremsstrahlung. The third and fourth terms are the results of 
the interference of the final state radiation from the quark and other radiation processes (initial 
state radiation and radiation from the gluon line). They contain both double hard processes and 
interferences between hard-soft and double hard processes in Fig. 2-(c) and (d). In the limit of soft 
gluon radiation, (1 — z) — > 0, these three terms can be neglected and we recover the result in the 
helicity amplitude approximation [13]. 



IV. MODIFIED FRAGMENTATION FUNCTION AND PARTON ENERGY LOSS 



Substituting Eq. (22) into Eq. (11), Eq. (5) and adding the gluon fragmentation processes, we 
have the semi-inclusive tensor from double quark-gluon scattering including the contribution beyond 
the helicity amplitude approximation, 

- 1 dz 

'Zh 



W D ' q 
dz h 



= ]T J 'dxHj${xp,q) j^D q ^ h (z h /z)^C A ^- 



dt\ 2ira s 



T q A g (x,x L ) + (1 - z)T«U(x,x L ) + g(l - z) 2 T^\x,x L ) 



+ (g — fragmentation) + (virtual corrections) . 



(23) 



where 



T qg (x,x L ) = I — 



dy^dy^il - e~ iXLP+y ^ )(1 - e - lXLP+{y )y{x+x L ) P + y - 



-(A\UQ)i + F+{yi)F + ^y^)Uy-)\A)0{-y2)°{y--yi) , 

T q f 1 \x,x L )= j ^-dy-dy^ ^ e - iXLP+ ^~ ~ y ^^ + e~ iXLP+y ^ - 2e~ iXLP+ ^~ - y ^ +y ^^ 



(24) 
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e i(x+ ^ )p+ ^i(A|^(0) 7 + ^( y2 -)^(yr)^(y-)l^>e(-?^)%- , (25) 



qy J 2n 



\(A\i> q (0) 7 + F a +(y;)F+°(y^Tp q (y-)\A)e(-y;)6(y- - y^) (26) 

are twist-four partem matrix elements of the nucleus. Apparently these parton matrix elements are 
not independent of each other. T q i g (x,x L ) has the complete four terms of soft-hard, double hard 
processes and their interferences. Therefore it contains essentially four independent parton matrix 
elements. T q ^ 1 \x,Xl) and T^ 2 \x,xl) are the results of the corrections beyond helicity amplitude 
approximation. But these two matrix elements are already contained in T^(x,xl)- 

During the collinear expansion, we have kept It finite and took the limit kx — * 0. As a conse- 
quence, the gluon field in one of the twist-four parton matrix elements in Eqs. (24)-(26) carries zero 
momentum in the soft-hard process. However, the gluon distribution xf g {x) at x = is not defined 
in QCD. As argued in Ref. [13], this is due to the omission of higher order terms in the collinear 
expansion. As a remedy to the problem, a subset of the higher-twist terms in the collinear expansion 
can be resummed to restore the phase factors such as exp(ixrP + y~) , where xt = (k T ) /2p + q~ z is 
related to the intrinsic transverse momentum of the initial partons. As a result, soft gluon fields in 
the parton matrix elements will carry a fractional momentum xt- 

Using the factorization approximation [13,14,17] we can relate the twist-four parton matrix ele- 
ments of the nucleus to the twist-two parton distributions of nucleons and the nucleus, 

T q A g (x,x L ) = ^-(l-e-^/^if^x + x^XTf^x^ + f^xXxL+x^f^XL+XT)] (27) 

XA 

where C is a constant, xa = 1/MRa, fq{x) is the quark distribution inside a nucleus, and f g N {x) is 
the gluon distribution inside a nuclcon. A Gaussian distribution in the light-cone coordinates was 
assumed for the nuclear distribution, p{y~) = /9oexp(?/~ /2i?^ ), where = \j2RaM / 'p + and M 
is the nucleon mass. We should emphasize that the twist-four matrix elements is proportional to 
1/xa = RaM, or the nuclear size [17]. 

Notice that the off-diagonal matrix elements that correspond to the interferences between hard- 
soft and double hard processes is suppressed by a factor of exp(— x 2 L /x 2 A ). This is because in the 
interferences between double-hard and hard-soft processes, there is actually momentum flow of 
xlP + between the two nucleons where the initial quark and gluon come from. Without strong long 
range two-nucleon correlation inside a nucleus, the amount of momentum flow x^p + should then 
be restricted to the amount allowed by the uncertainty principle, 1/i?^ ~ p + /RaM. Similarly, the 
other two parton matrix elements in Eqs. (25) and (26) can be approximated as 

T^\x,x L ) = ^{[^{x + x l )xt^{xt) + ^{x){x l +xt)^{x l +xt)} e~ x ^ 

- 2f q A (x)(x L + x T )fg{x L + x T )} , (28) 
T£*\x,x L ) = —tf{x){x L + x T )f^{x L + x T ) . (29) 

XA 

From the above estimate of the matrix elements, both T q ^ j (x, Xl) and T q ^ 1 \x, xl) contain a factor 
1 — e~ x i/ x A because of the LPM interference effect. Such an interference factor will effectively 



9 



cut off the integration over the transverse momentum at xl ~ x a in Eq. (23). As we will show 
later in the calculation of the effective energy loss, the integration with such a restriction in the 
transverse momentum due to LPM interference effect will give rise to a factor 1/xa in addition to 
the coefficient f A (x) / xa ■ Consequently, contributions from double scattering in Eq. (23) that are 
associated with T a (x,xl) and Tqg 1 ^ (x, x£) will be proportional to R A f A (x). These are the leading 
double scattering contributions in the limit of large nuclear size. On the other hand, the third term 
T^ 2 \x,xl) in Eq. (23), which does not contain any interference effect, will only contribute to a 
correction that is proportional to RAf A (x). In the limit of a large nucleus, A 1 / 3 » 1, we will neglect 
this term in the double scattering processes. It is interesting to point out, however, that the physical 
process associated with this term is totally responsible for the non- vanishing photon spectra in QED 
bremsstrahlung which otherwise vanishes in the helicity amplitude approximation. As we can see, 
the leading correction beyond helicity amplitude approximation comes from the interference between 
this process and other radiation processes that contribute to the leading result in the first term. 

The virtual correction in Eq. (23) can be obtained via unitarity requirement similarly as in 
Ref. [13]. Including these virtual corrections and the single scattering contribution, we can rewrite 
the semi-inclusive tensor in terms of a modified fragmentation function D q ^h(zh, M 2 )i 



dW, 



dz h 



]T J dxf q A (x lf i 2 I )Hj°J(x,p,q)D q ^ h (z h ,n 2 



(30) 



where f A (x, fjij) is the quark distribution functions which in principle should also include the higher- 
twist contribution [18] of the initial state scattering. The modified effective quark fragmentation 
function is defined as 



D^ h {z h ,f) ee D^ h {z h ,f) + jf |f g f f [A 7g _ 



j(z, x, x L jV)D q ^ h {zh/z) 



+ Aj q ^ gq (z 7 x,x L ,lT)D g ^ h (zh/z)] , 



(31) 



where D q ^h{zh, V?) and D g ^h{zh, M 2 ) are the leading-twist fragmentation functions. The modified 
splitting functions are given as 

1 + z 2 



^Iq^qg (z, X, XL , £j<) 



.(!"*) + 

A7 M ,(z,i,i l ,4) = A~f q ^ qg (l - z,x,x L J T ) 

1 



< {m) {x, x L ) + 6(1 z)AT A ^(x,£ 2 T ) 



2ixu s Ca 



ATf™)(x,4) 



dz- 



1 - 



27^(^)1^ - (1 + z 2 )T A } m \x,x L ) 
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T$ m \x,x L ) ee T A (x,x L ) + (1 - z)T A ^(x,x L ) . 



(32) 

(33) 
(34) 
(35) 



The above modified fragmentation function is almost the same as in the previous calculation 
with helicity amplitude approximation, except that the twist-four parton matrix element T a (x,xl) 
is replaced by a modified one T qg m \x,XL) in Eq. (35). One can then calculate numerically the 
modified fragmentation function as in Refs. [12,13]. To further simplify the calculation, we assume 
xt <C xl <C x. The modified parton matrix elements can be approximated by 



C 



T A ^(x,x L ) « —(1 - e- x ^)f A (x) 



xa 



(36) 
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where C = 2CxtJ^{xt) is a coefficient which should in principle depends on Q 2 and xt- Here we 
will simply take it as a constant. The new correction term in this calculation is thus negative in 
the modified splitting function. This will reduce the nuclear suppression of hadron spectra at large 
values of z and thus reduce the effective quark energy loss. 

Because of momentum conservation, the fractional momentum in a nuclcon is limited to xl < 1. 
Though the Fermi motion effect in a nucleus can allow xl > 1, the parton distribution in this 
region is still significant suppressed. It therefore provides a natural cut-off for in the integration 
over z and It in Eq. (31). Shown in Fig. 3 are the calculated nuclear modification factor for the 
quark fragmentation function Da(z, Q 2 )/Djy(z, Q 2 ) inside a nucleus with A = 100. In this numerical 
evaluation, we have taken C — 0.006 GeV 2 which was fitted to the HERMES experimental data [12]. 
The dashed curve is for the modified fragmentation function in the hclicity amplitude approximation 
and the solid curve is obtained with the new correction term. Apparently, the new correction term 
reduces the nuclear modification, though the reduction is not very significant. 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



FIG. 3. Calculated nuclear modification factor for the quark fragmentation in a nucleus (A=100). The 
solid line is the current calculation with the new correction term. The dashed line is the previous result with 
helicity amplitude approximation. 



Similarly, we can also calculate the effective quark energy loss, which is defined as the energy 
carried away by the radiated gluon, 



rv 2 d( 2 f 1 a s 
(Azg}(x B ,^ 2 ) = J J dz^-zAj q ^ gq (z,x B ,x L J T ) 



C A a 2 8 



Jo "t Jo 



N c Jo 4 ^0 ^ ' ' f£{x B ,tf) 

We separate the parton energy loss as two parts 

(Az g )(x B ,H 2 ) = (Az g ) he li(x B ^ 2 ) + (AZg) cor r(xB,H 2 ) 



(37) 



(38) 



where (Az g )i ie ii{xB, M 2 ) is the leading quark energy loss with helicity amplitude approximation [13], 
and {Az g ) corr (x Bl M 2 ) is the new correction to the quark energy loss in this calculation. Using the 
approximation for the modified twist-four parton matrix elements in Eq. (36), we have 
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(41) 

where x M = [i 2 /2p + q~ z(l — z) = Xb/z(1 — z) if we choose the factorization scale as fi 2 = Q 2 . When 
xa <C xb <C 1 we can estimate the leading quark energy loss roughly as 



(Az g ) hell {x Bl f) « C ° Aa " X 2 B 2V^l3\n 1 —^~ 1] , . 12) 

iV c IcJ X %B 



(Az g ) corr (x B , M 2 ) « Q5fi ~ ^ + o] • ( 43 ) 



Since xa = 1/MRa, both of the energy loss (Az g )heli and (Az g ) corr depend quadratically on the 
nuclear size. Adding them together, we have 



Ca 2 x B ^ r _. l-2x B 5 
A?TQM^ [51n ^ 2 ] 



(Az g )(x B ^ 2 ) « -^T^F^ln— — - ^ - U) 

5 



- -{Az 9 ) keli (i B , m ), (au < < 1) . (45) 

The new correction thus reduces the effective quark energy loss by approximately a factor of 5/6 
from the result with helicity amplitude approximation. 

V. SUMMARY 

We have extended an earlier study [13] on gluon radiation induced by multiple parton scattering 
in DIS off a nuclear target with a complete calculation beyond the helicity amplitude (or soft radi- 
ation) approximation. Working within the framework of the generalized factorization of twist-four 
processes, we obtained a new correction to the modified parton fragmentation functions. Such a new 
correction essentially results in a new term in the modified splitting function which is proportional 
to (1 — z). In the limit of helicity amplitude approximation (1 — z) — * 0, this term vanishes and we 
recover the early results [13]. 

The new correction we obtained in this paper comes from the gluon radiation process (residual 
final state radiation from a quark after incomplete cancellation by the initial state radiation) that is 
actually responsible for the photon radiation in QED. However, the leading correction beyond the 
helicity amplitude approximation does not come from this process itself. Rather, it comes from the 
interference between this process and the other gluon radiation processes that are responsible for 
the result in the helicity amplitude approximation. Though it is not dominant for induced gluon 
radiation in QCD, it still make a finite contribution to the modified fragmentation function for a 
quark propagating inside a nuclear medium and to the effective quark energy loss. We found that 
it reduces the effective quark energy loss by a factor of 5/6. 
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APPENDIX 

In this Appendix we will list our complete calculation of quark-gluon double scattering in detail. 
There are total 23 cut diagrams which are illustrated in Figs. 4-14. 



xpi 






Ap V i i i 



Ap 



FIG. 4. 



In Fig. 4 there are three possible cuts(central cut, left cut and right cut). We get 

di% a s „ 1 + z 2 



H A P ~AV .2/1 .3/2 i k T,x,p,q,z) = J -JT^ C F~^ 



T 
2-KOi. 



where 



x -jj— lA P -i(y .3/i ,3/2 ,^T,k T , x, p,q,z), 



(46) 



lA P -i,c(y~,yi,y2^T,k T ,x,p,q,z)=e^ +x ^ » +«°" ^ -*> >9(-y^)9(y- - yf) 

x (1 - e - iXLP+v *)(l - e ~ lXLP+(v ~- y i*>) . (47) 

i Ap -i,L(y'>yi>y2>iT > k T ,x,p ) q,z) = -e^+^P+^+^+fcr-^)^- _ y - )0 ^- _ yT) 

x ( 1 _ e -ixL P +(y--y-)^ ( 48 ) 

I Ap . hR (y-,y^,y^,l T ,k T ,x,p,q,z) = ~e^ +x ^ + ^ + ^ + ^-y^9(-y^)9(y^ - yf) 

x (1 - e -i!EtP+! '0 . (49) 
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FIG. 5. 

In Fig. 5, there are two different cuts, central or left. So we obtain, 



H A P -2(y ,Ui ,V2 ,kT,x,p,q,z) = j dC T 



2 £ T ■ (It - (1 - z)kr) a 

T~ 

27ra s — 



e T {i T - {i - z)k T ) 2 27r 



... r , C A 1 + z 2 



2 ' 1 - z 



x -j-j—lA P -2(y ,yi,y 2 ^T,kr,x,p,q,z), 



_ „i(x+a;z,)p+3/ +ia:DP + (j/ 1 -J/ 2 W_ 



2/2 )% - ih ) 



(50) 
(51) 



lA P -2,c(y ,2/i ,2/ 2 ,^T,k T ,x,p,q,z) = e % 

x [e 

/A P -2,L(y-, 2/f, ife ,4r, *r, ar,p, 9, z) = e ^ + ^^ + ^ + ^ + ^-^)%- - yf)%r ~ y 2 ") 
where x° D = kx/2p + q^ . 



xpi 




> v » 

x 3 p+k T 


— *■ V- » 

]g x 2 pH-k T "x 











Ap 



FIG. 6. 

As for the central cut and right cut of Fig. 6, we obtain 

l T ■ {t T - (1 - z)k T ) a 



H A P -3(y >Vi ,2/2 ,kr,x,p,q,z) = / d£ T "^ ^ — ^7^7 t£ [C f ~ ^) 



CV 1 + Z 2 



4(£ T - (1 - «)*r) 



2 2vr 



2 ' 1 - z 



2ira s 



X -j-j—lA P -a{y ,Vi ,y 2 Jr,kT,x,p,q,z) , 



(53) 



/Ap-3,c(2/^2/r,2/ 2 ^T,fcT,x,p, 9 ^ ^ 

x [ e - ia; £P + 2^ _ e — «:iP + (3r-3/r+!';r)l j (54) 

7^-a 1 K(l/-,tf 1 -,l/ 2 -,ir,*r,a:,P,?,*) = e^ +x ^ + ^ + ^ + ^^9{-y^)9{y^ - y^) 

x [ e - ,; ( a;f i- 2: D)p + (ar"-y2")- !;;i: iP + yr _ e - ix LP + v2^ _ ^55) 
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FIG. 7. 



There is only one cut (left cut) in Fig. 7 with the contribution, 



Tf D i - - - ; ^ [ di T a 

H A P -i(y ,yi,V2 ,k T ,x,p,q,z) = j ^--^ 



\ + z z 
1 - z 



Ia p -a,l{.v ,Vi ,y 2 ,lT,kr,x,p,q,z 



x ~^- ll Ap-i(y ,Vi ,V 2 ,ir,k T ,x,p,q,z) , 

= _ e i(.x+x L )p+y-+ix D p+(y--y-) e ^ y - _ y -)Q( y - 
x e i i x °D- x D)p + (Vx -Vz) e ~ix L p+ (y- -y~) ^ 



xp 



Ap 



x 3 p+k T 



x 2 p+k T 3| -xp 



Ap 



As for Fig. 8, we get 



FIG. 8. 



-rrD 



H Ap-b(v >Vi .2/2 ,k T ,x,p,q,z) 



a s 1 + z z 
— C F - 

1 — z 



{l T - (1 - z)k T ) 2 2vr 
x ^j^lAp-5 (y~ , Vi , y 2 ,?t, k T , x, p, q, z) . 
I Ap - 5 ,c(y-,yi,y 2 jT,k T ,x, P ,q,z) = e^ + ^ + v- +ix ™ + ^-^e{-y 2 )9{y- - 

x e ~ixLP + {y~ ~Vi) e ~ix L p + y~ 
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1 . . & 






X 3P +k To 


|x 2 p+k T 


"X 











Ap 



FIG. 9. 



The contribution from Fig. 9 is 

H A P ~6(y >Vi iV2 ,k T ,x,p,q,z) = J -p- ^ C F — 



T 

Ap—6 iy~ j V\ j V2 ; ^-T-, k T , x,p, q, z) , 



(60) 




FIG. 10. 

As for the processes in Fig. 10 we have three possible cuts. Thus the contributions can be written 



as 



Ap-7 



civ >Vl >2/2 ,k T ,x,p,q,z) 



{It - hr) 



2 2^ Ca l - 7 



2-7ra, 



x —j^—lA P -7,c(y ,Vi ,y 2 ,tT,kr,x,p,q,z), 

ix D p+y~ /(1-z) _ e -ixLP + V2~l 

n+C / -_«-1/('i_»1 _iT,^+^r. 



x e" 



H Ap-7,L(R)iv >Vl >V2 > k T,x,p,q,z) = J ^^ G Aj 



x [g^DP (W -» x _ g-«B£P + (l/ "Si )l 

2 - 1 + z 2 



(62) 



(63) 



27TQI 

x -j^-lA P -7MR)iy~ 'Vi ^2 JT,k T ,x,p,q, z), (64) 
7^-r,£(»-,y 1 - > » a -^r,fcr,i,p,g,z) = -e^+^^+^^r-^)^- _ y -)e(y- - y ~) 
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i(l-z/(l-z))x D p + (y 1 -y 2 ) 



X e 

n + 



x [1 - e- lXLP (y -«i )], 




FIG. 11. 



There are three cuts in Fig. 11 and we have 



H A P -s(y >Vl 1 2/2 , k T ,x,p,q, z) 



2 l T ■ {l T - k T ) a s C A l + z 2 



l 2 T {l T -k T ) 2 2tt 2 1-z 



x lAp-${y ,y 1 ,y 2 ,£ T ,kT,x,p,q,z), 
lA P - 8 ,c(y-, Vi, Vi^T, kr, x,p, q, z) = -e^ + *^ + y- +i *°P + ^-y^0{-yz)9{y--yi 

x J e «DP + J/^/(l-z) _ e -ix L p + y~-i 

x [i - e -^p + (^^D] 5 

x j e -^i,p + (y _ -ar) _ e i2; DP + (s/ _ -yr)/( 1 ^ z )] j 
x [ e _i:Cj: 'P + J'r _ e «DP + s 2 "/( 1 - z )] 




FIG. 12. 
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The contributions from the three cuts in Fig. 12 can be read as 

d _ _ f in2 It ■ (g - k T ) a s C A l + z 2 

H Ap _ 9 (y, Vl ,y 2 ,kT,x, P , q ,z) = Jd£ T ^ _ ^ — 

x -^7A P -9(y" ,Vi \V% ,tT,k T ,x,p,q,z) , 

I Ap ^ iC (y-,y^,y 2 ,£ T ,k T ,x,p,q,z) = _ e K*+^)P + iT+Wt B r- ^e(-y^)6(y- - yf) 

x dp + ~(y~ -vl ") / '(1- z) _ g-»«-E,P + (v~-vr)] 

x [1 - e- tol3,+ ^] , 

T^- 9 ,L(»-,yr,y2^T,fcr,a:,p,«,«) = -e^^"^^^-^)^- - y r)%r - 

x { e ~ iXLP+ ( y ~~Vi) — e ix £>p + (i/ _ -3/i~)/(l-2)i 

x r e -ixtp+K~ _ gia;£5p+jf 2 "/(l-z)j 



Ap 



q 






3 , f 


P'' o 


t 






"X 















Ap 



FIG. 13. 



In Fig. 13 there are two possible cuts (central or left). We have 



-rrD 



H A P -io,c(y iVi ,2/2 , k T,x,p,q,z) = / dt 



(It - kr) ■ (It - (1 - z)kr) a s C A 1 + z 1 
' (l T - kTY(€r - (1 - z)k T ) 



2 2tt 2 1 - z 



x -jj—lA P -lO,c(y ,Vi ,2/ 2 ,^T,k T , x, p, q, z) , 



Ia p -W,c(V >Vl ,2/2 jT,k T ,x,p,q,z) 



_ i(x+x L )p + y +ix D p + (y 1 -y 2 ) nr_ 



2/2 )% - 2/i ) 



x e ~ixLP y 2 \ e ix DP + {y -y ± )/{l-z) _ & -ix L p + (y -y 1 )j 



— D 



H Ap-io,r(v ,2/i ,2/ 2 , k T,x,p,q,z) = / <2£ 



£ T • (4r - «M a s CU 1 + z 2 



27ra, 



x -jj—lAp-ioMv >Vi >V2 JT,k T ,x,p,q,z) , 
x [e _i ( :Cr>_a: D)p + (fi~-3'2 )-»-Lp + i'r' 



-i(l-z/(l-z))a; D p + (i/ 1 -y 2 )-ix L p + y 2 
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FIG. 14. 



In Fig. 14 we can make the central cut or the left cut and we obtain 

fl2 (g. - k T ) ■ (lr - (1 - z)k T ) as_ C A l + z 2 
' (f T - k T ) 2 {f T - (1 - z)k T ) 2 27r 2 1 - z 



-rrD 



H A P -n,c(y iVi ,V2 ,k T ,x,p,q,z) = / dt 



x '—^j— 1 A P -n,c(y ,Vi ,V 2 Jr,k T ,x,p,q,z) , 

i Ap -n,c{y~,yi,y2JT, **, x, P , q , z) = e^+^^+^+ci/r-^-) e(- y -)e( y - - y -) 

x e -i*LP + {y~ -y±) ^ e ix D p + Vz /{l-z) — e - ix LP + Vz^ 



-ttD 



H A P -u,R(y '2/2 , k T,x,p,q,z) = / 



£ T • (£ T - zkr) a s C A 1 + z 2 



4(£ T -2fc T ) 2 2tt 2 1-2 



x ^^^ p -ii,l(2/ ,y 2 ,^T,k T ,x,p,q,z) , 

-i{x d -x a D )p + {y~ -y~ )-ix L p + (y~ -y~ ) 

n~^~ ( ii — ii \ — if r ( 



x e 



-i(l-z/(l-z))x D p + (y 1 -y 2 )-ix L p+(y -y x )j 
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